Draft version March 8, 2013 

Preprint typeset using lATgX style emulateapj v. 11/12/01 



o 



Oh 



AGNS AS MAIN CONTRIBUTORS TO THE UV IONIZING EMISSIVITY AT HIGH REDSHIFTS: 
PREDICTIONS FROM A A-CDM MODEL WITH LINKED AGN/GALAXY EVOLUTION 

E. GiALLONGO, N. Menci, F. Fiore, M. Castellano, a. Fontana, A. Grazian, L. Pentericci 

INAF - Osservatorio Astronomico di Roma, via di Frascati 33, 1-00040 Monteporzio, Italy 

Draft version March 8, 2013 

ABSTRACT 

We have evaluated the contribution of the AGN population to the ionization history of the Universe 
based on a semi- analytic model of galaxy formation and evolution in the CDM cosmological scenario. 
The model connects the growth of black holes and of the ensuing AGN activity to galaxy interactions. In 
the model we have included a self consistent physical description of the escape of ionizing UV photons; 
this is based on the blast- wave model for the AGN feedback we developed in a previous paper to explain 
the distribution of hydrogen column densities in AGNs of various redshifts and luminosities, due to 
absorption by the host galaxy gas. The model predicts UV luminosity functions for AGNs which are in 
good agreement with those derived from the observations especially at low and intermediate redshifts 
{z ~ 3). At higher redshifts {z > 5) the model tends to overestimate the data at faint luminosities. 
Critical biases both in the data and in the model are discussed to explain such apparent discrepancies. 
The predicted hydrogen photoionization rate as a function of redshift is found to be consistent with that 
derived from the observations. All that suggests to reconsider the role of the AGNs as the main driver 
of the ionization history of the Universe. 

Subject headings: galaxies: active — galaxies: formation — galaxies: evolution 



1. INTRODUCTION 



Q The assessment of the thermal history of the intergalac- 
^ ic medium at high redshifts is fundamental to under- 
C/5tand the physical processes involved in galaxy formation 
C^nd evolution, including triggering and flueling of Active 
' Galactic Nuclei (AGNs). 

The ionization state of the intergalactic medium (IGM) 
^s a function of the hydrogen ionizing UV background 
0(UVB) conceivably produced by sources like high redshift 
l/^alaxies and AGNs. Star forming galaxies are more com- 
OVion in the high redshift universe, and in principle can 
C^e responsible for the reionization of the IGM starting at 
> 7. However, even a first order estimate of their con- 
Qk ribution to the ionizing UVB has two main sources of 
C^ncertainty. The first is connected with the increasing dif- 
^^culty in evaluating the faint end slope of the UV luminos- 



have been made to derive the escape fraction of UV ioniz- 
ing photons both at low redshifts from space and at high 
redshifts z = 2 — 4 from ground based observations. How- 
ever, there is little evidence supporting a scenario where 
enough ionizing photons escape from galaxies at interme- 
diate redshifts, and most attempts gave only upper limits 
on the escape fractions. 



Indeed rece nt measurements by |Cowie, Barger, 
Trouille (2009|) provided a stringent 2a upper limit of 0.^ 



^ty function especially at the hig hest redshifts ( Bouwens et 
.r4l^ ||2011J|G razian et al. ||2011|). The second is even less 
^<Jinderstood and is connected with the poor knowledge of 
J-ihe average escape fraction of ionizing Lyman continuum 
photons from the interstellar medium of each galaxy. On 
the other hand, the AGN contribution is mainly affected 
by the poor knowledge of the faint end slope of their lu- 
minosity function especially at high redshifts due the the 
lack of deep AGN surveys at various wavelengths (see e.g. 
Shankar & Mathur (2QQ7[)). The observed escape fraction 



for the relative escape fraction of galaxies at0.9<z<1.4 
suggesting that galaxies are not the dominant contribu- 
tors to the ionizing UV background at 2; < 2. At high 
redshifts the situation is more difficult from the observa- 
tional point of view. Direct measures of the Lyman con- 
tinuum flux from galaxies at z > 4 are increasingly diffi- 
cult because of the sharp increase in the number density 
of intervening Lyman Limit absorptions due to optically 
thick clouds in the IGM. For this reason the search has 
been focused to the range z ~ 3 — 4, but there a discrep- 
ancy is found by different teams using different methods. 
Steidel and collaborators support the detection of appre- 
ciable Lyman limit flux es caping at least from the UV 
brightest galaxies at z ~ 3 (Steidel, Pettini, & Adelberger 
] |2QQ1 ). More recent analyses suggest average es cape frac- 
tions of order 10 - 15% (Shapley et al. ||2QQ6[ [Nestor et 



al. 2011), although t heir sample can be contaminated 



of ionizing photons in the majority of the bright AGNs is 
high and of the order of unity at the Lyman continuum, 
although the value for the faintest AGNs is unknown. 

Since the apparent number density of bright QSOs and 
AGNs is rapidly decreasing outside the redshift interval 
2 < z < 3 it is usually assumed that star forming galaxies 
can give a considerable contribution to the ionizing flux 
which could become dominant at z < 1 and at z > 3 
^Madau ||1 99l! 'Giallon go, Fontana, Madau ||1997[ [Haardt 
L Madau i.l996, 2012, ). In this context several attempts 



by foregrou nd objects (Vanzella et al. 2010 Boutsia et 



al. ||2011 ). Other teams using both spectroscopic and 
very deep broad or narrow band imaging gave only upper 
limits in the r ange < 5 — 15% (Giallongo et al. 2002] 
Fernandez- So to, Lanzetta, fc Chen ||2003| [Vanzella et al.| 
m2010b^ Boutsia et al. 2011). To summarize, the recent 



trends appear to favor low escape fractions from relatively 
bright galaxies, where the recent limits are below 1% at 
low redshifts and below 10% at z ~ 3 — 4, and cast serious 
doubts on any redshift evolution in the galaxy ionizing es- 
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cape fraction. This leaves open the possibihty that AGNs 
could still play a leading role in contributing to the cos- 
mological ionizing background. 

The traditional view based on relatively bright QSO sur- 
veys in optical and X-ray bands lead to the standard view 
of an evolutionary ionizing QSO background increasing 
from the local value out to 2: = 2 — 3 and then quick ly 



decreasing (see e.g. Madau, Haardt, & Rees ([1999)). 



Previous estimates were based on simple parametric ex- 
trapolations both in luminosity and redshifts of the AGN 
luminosity function without any link to sp ecific physical 



processes (Madau, Haardt, & Rees 
cent deep otpical QSO surv eys aUtT 



2008 



1999[ ) . However re- 
= 3 (Siana et al. 



Fontanot et al. |[2QQ7| ) and z = 4 ( [Glikman et al. 



2Qlip are showing the presence of a considerable number of 
previously unknown faint AGNs producing a rather steep 
luminosity function. The inclusion of X-ray detection in 
the selection methods enables to extend the k nowledge of 
the lu minosity function to even fainter limits (Fiore et al. 
] |2Q12[ ) of the order of L(2-10 KeV) > 10^^ erg s"^ out to 
z > 4. 

The presence of this faint population is changing our es- 
timate of the AGN contribution to the ionizing UV back- 
ground although the selection of faint AGNs at the high- 
est redshifts z ^ 6 becomes difficult with the current in- 
strumentation. Given the current observational limits, an 
insight into the origin of the ionizing flux at high red- 
shifts can be gained from theoretical modelling. The latter 
should account for the formation and evolution of galax- 
ies in a cosmological context, and for the co-evolution of 
the AGN population. Given the complex sub-grid physics 
involved in the AGN activity over a wide range of cosmic 
times and of galactic masses, semi-analytic models (SAMs) 
constitute a powerful tool to provide a statistically rele- 
vant sample of simulated galaxies and AGNs. Such mod- 
els connect the physical processes involving gas and star 
formation to the merging histories of DM haloes collapsed 
from the primordial density field. While several SAMs in- 
cluding the growth of supermassive Black Holes have been 
presented i n the literature (starting from Kauffmann fc 
Haenhelt | (|2QQQD; see alsolCroton et al. | PQQ6D; |Bower " 



et al I ( |2QQ6P ; [Marulh et al. | ( |2QQ8p ; |Menci et al.' | ( 120061 
2008p) , addressing the problem of the ionizing background 
requires a physical description of the escape fraction in 
galaxies. Such a problem is deeply connected with that 
concerning the feedback of the AGNs onto the interstellar 
gas. In fact, since even a tiny amount (column densities 
Nh ~ 10^^ cm~^) of galactic gas would suppress the ra- 
diation at the Lyman limit escaping from the galaxy, esti- 
mates of the escape fraction translate into computing the 
fraction of photons that are emitted along directions where 
the interstellar gas of the g alaxy h as been depleted. In our 
previous work ( Menci et al. |2008|) , we have developed our 
SAM model (Rome-SAM, hereafter R-SAM)) to describe 
the expansion of the blast wave produced in the galaxy 
interstellar medium by AGN-driven outflows. Here we use 
the same model to derive the escape fraction of photons in 
active galaxies as the fraction of photons emitted along di- 
rections where the AGN-induced blast wave has expelled 
the galactic gas. When coupled with the integrated lu- 
minosity function of AGN self-consistently derived in the 
same model, this allows us to compute the ionizing back- 



ground produced by AGNs within a cosmological model of 
galaxy formation. We stress that the basic quantities pre- 
dicted by our computation have been extensively tested 
against observations in our previous works. On the one 
hand, the predictions concerning the growth of supermas- 
sive black holes and the corresponding evolution of the 
luminosity function have been compared against a wide 
range of different observations, from optical to X-rays, for 
SMBHs masses and AGN luminosities in a wide range of 
cosmic times ex t ending up to z=6 f or the brightest QSOs 
dMenci et al | ( |2003[ |2Q04[ |2006j |2008[ ); see also [Fiore 
|et al. I (|2012p); on the other hand, the blast-wave model 
for the AGN feedback allowed us to derive -within a full 
cosmological context and in connection with the galaxy 
properties- the distributions of column densities in active 
galaxies, as a function of redshift and of the AGN lumi- 
nosity, which we successfully tested against observations 
(Menci et al. 2008). In this paper we extend the analysis 
of the UV AGN luminosity function to faint objects and 
adopt the blast- wave mechanism to predict the amount of 
UV photons escaping from the AGN host galaxies. 

The paper is organized as follows: in Sect. 2.1 we re- 
call the basic properties of our semi-analytic model for the 
evolution of galaxies and AGNs. In Sect. 2.2 we describe 
how our blast-wave model for the AGN feedback (devel- 
oped in previous papers) allows to compute the luminosity- 
dependent escape fraction of ionizing photons emitted by 
AGN. The evolution of the luminosity function of AGNs 
predicted by our model is compared with data in Sect. 
3; the comparison with the current dataset is discussed 
taking into consideration possible biases in the model and 
in the data. In Sect. 4 we convolve the predicted in- 
trinsic luminosity function of the AGN with the escape 
fraction derived in our model to derive the emissivity and 
the photoionization rate predicted by the model at various 
redshifts, and compare the latter with existing data; the 
contribution of AGNs with different luminosities to such 
observables is also shown and discussed. The final Sect. 5 
is devoted to a summary and to conclusions. 

2. AGNS EMISSION IN HIERARCHICAL GALAXY 
EVOLUTION 

The semi-analytic model R-SAM we develop and use 
connects the AGN evolution with the evolution of gal axies 
in a cosmological framework (see Menci et al. (2008) and 
references therein). 

2.1. Hierarchical Galaxy /AGN Evolution in the R-SAM 

Galaxy formation and evolution is driven by the col- 
lapse and growth of dark matter (DM) haloes, which orig- 
inate by gravitational instability of overdense regions in 
the primordial DM density field. This is taken to be a 
random, Gaussian density field with Cold Dark Matter 
(CDM) powe r spectrum expected wi thin the "concordance 



cosmology" (Spergel et al. 2007) for which we adopt 
round parameters Ua = 0.7, = 0.3, baryonic density 
= 0.04 and Hubble constant h = 0.7 in units of 100 
km/s/Mpc. The normalization of the spectrum is taken to 
be (78 = 0.9 in terms of the variance of the field smoothed 
over regions of 8 Mpc. The merging histories of DM 
haloes and sub-haloes is followed through a Monte Carlo 
simulation of the collapse and subsequent merging history 



Giallongo et al. 
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of the peaks of the primordial density field, which enables 
us to generate a synthetic catalogue of model galaxies and 
of their past merging history. The physics of baryons (gas 
radiative cooling, disk formation, quiescent and impulsive 
star formation, stellar and AGN feedback) is connected to 
the merging histories of D M haloes as describe d in detail in 
our previous papers (e.g., [Menci et al. ( [2QQ5| )). Note that 
the impulsive star formation implemented in the model is 
triggered not only by (major and minor) merging events 
but also by fiy-by events between galaxies. The R-SAM 
also includes a treatment of SMBHs growing at the cen- 
tre of galaxies by interaction-triggered infiow of cold gas, 
following the physical model propos ed byjCavaliere fc Vit 



torini 



2QQQD and implemented in [Menci et al. j ( |2QQ3| 



2006, 2008p . The accretion of cold gas is triggered by 



galaxy encounters (both of fly- by and of merging kind), 
which destabilize part of the galactic cold gas mass rUc 
by inducing loss af angular momentum. The fraction of 
cold gas accreted by the BH in an interaction event is 
computed in terms of the variation Aj of the sp eciflc an- 
gula r momentum j ^ Gm/vd of the gas, to read (Menci et 
al. [[2003) face ^ 10~^ |Aj/j|. For minor merging events 



and Tor the encounters among galaxies with very unequal 
mass ratios m'/m <C 1, that dominate the statistics in 
all hierarchical models of galaxy formation, the accreted 
fraction takes on values 10~^ < face ^ 10~^. The average 
amount of cold gas accreted during an accretion episode 
is thus Arriacc = facc'^ci and the duration of an accre- 
tion episode, which provides the timescale for the QSO or 
AGN to shine, is assumed to be the crossing time r of the 
destabilized cold gas component. The time- averaged bolo- 
metric luminosity so produced by a QSO hosted in a given 
galaxy is then provided by 

^ r]c^ Arriacc 



We adopt an energ y- conversion efliciency 7^ = 0.1 (see[Yu| 
fc Tremaine ( 2002[ )), and derive luminosities in the various 
ban ds on adopting standa rd spectral energy distributions 
as in [Marconi et al. ( [2004 ) and Marconi et al. (in prepara- 
tion )^nie"^MBTrniass rriBH grows mainly through accre- 
tion episodes as described above, besides coalescence with 
other SMBHs during galaxy merging. As i nitial condition, 
we ass ume small seed BHs of mass 10^ Mq (Madau fc Rees 
|2001 |) to be initially present in all galaxy progenitors; our 
results are insensitive to the speciflc value as long as it is 
smaller than some 10^ Mq. Due to the crucial role of AGN 
feedback in affecting the star formation history of massive 
galaxies, this part of the model has been extensive tested. 
The predicted QSO luminosity function and X-ray lumi- 
nosity function and emissivity have been checked against 
observation over a wide range of redshift < z < 5 and 
bolometric luminosities 10^^ < L < 10^^ erg s~^ ( Menci et 
The consistency of the modeling of the 



1120031 

Nra 



2004) 



back with the colour distribution of galaxies, in- 



cluding a descripti on of the bimo dal appearance, has been 

(20061). 



discussed in ^Menci et al. 



2.2. Blast-Wave Model for AGN feedback and the escape 
fraction of UV ionizing radiation 

Fast winds with velocities up to 10~^c are observed in 
the central regions of AGNs. They likely originate from 
the acceleration of disk outflows by the AGN radiation 



fleld (see Begelman] ( 2003 ) for a review) , and affect the 
environment in the host galaxy and beyond, leaving im- 
prints out to large scales of some 10^ kpc in the intracluster 
medium (ICM). Clear examples of outflows in local AGNs 
are those observed on galactic scales in the infrared and 
optical bands in Mrk 231 a nd Mrk 573 ([Feruglio et al. 
2Q10| [Fischer^^etjl~[ [2010| [Rupke fc Veilleux [[20 1 1|. A 



detailed model for the transport of energy from the inner. 



outflo w region to larger scale s ha s been developed by[Cav- 
alier e^Lapi, fc Menci (2002) and Lapi, Cavaliere, fc Menci 
y ( |2005| ). Central, highly supersonic outflows compress the 
gas into a blast wave terminated by a leading shock front, 
which moves outwards with a lower but still supersonic 
speed and sweeps out the surrounding medium. The key 
quantity determining all shock properties is the total en- 
ergy AE injected by AGNs into the surrounding gas. This 
is computed as 



= eAGNr]c^ Arriacc ^agnLt 



(2) 



for each SMBH accretion episode in our Monte Carlo 
simulation; the value of the energy feedback efliciency for 
coupling with the surrounding gas is taken as cagn = 5 x 
10~^, which is consistent with the values required to match 
the X-ray properties of the ICM i n clusters of galaxies (see 
Cavaliere, Lapi, fc Menci (20021)). This is also consistent 
with the observations of wind speeds up to ^ 0.1c in the 
central regions, which yield cagn ^ Vw/2c 0.05 by mo- 
mentum conservation between photons and p articles (see 



Chartas et al. (2002[); Pounds et al. (2003[)); this value 



has also been adopted in a number of simulations (e.g., Di 
Matteo, Springel, fc Hernquist (2005) ) and in our R-SAM 



(.e.g. 



Menci et al. 



The blast expani 



[ po6p j: 

ids into the ISM or the ICM as described 
by hydrodynamical equations that include the effects not 
only of an initial density gradient, but also those of up- 
stream pressure and DM gravity , clearly important quan- 
tities in galaxies, as discussed by [Lapi, Cavaliere, fc Menci 
Ij ( |2005| ). The solutions show in detail how the perturbed 
gas is conflned to an expanding shell bounded by an outer 
shock at radius Rs{t), which sweeps out the gas surround- 
ing the AGN. An analytic expression for Rs{t) has been 
derived by the above authors in the case of shock expan- 
sion in a gas with power-law density proflle p ex r~^; in 
the relevant central regions of the galaxies (within 1/2 the 
scale-length r^ of the ex ponential disk) the e xpansion of 
the shock is given by (see [Menci et al. [ ( [2008| )) 

Rs{t)^rdMt/td (3) 

where A4 is the Mach number. This is related to the en- 
ergy AE (eq. 2) injected by the AGN (relative to the 
initial thermal energy content E oc rrir. of the g:alactic ^as) 
by the relation = l-\-AE/E (Lapi, Cavaliere, fc Menci 
2005). Thus the expansion velocity of the bast wave is 
directly related to the AGN luminosity. 

The above blast-wave model for the AGN feedback al- 
lows to self-consistently compute the escape fraction of 
ionizing photons. In fact, the typical galactic column den- 
sities are in the range Nh ^ 10^^ ^ 10^^ cm~^ and effec- 
tively suppress the escape of ionizing radiation along the 
line-of-sights where the galactic gas has not been swept 
out by the blast wave produced by the AGN. Thus, the 
escape fraction can be computed as f{L) = 0^/2 tt, where 
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uj is the solid angle subtending the intercept of the expand- 
ing blast wave with the galactic disk. This is illustrated in 
fig. 1, which shows how only ionizing photons propagat- 
ing along the directions (indicated with £i in the figure) 
within the solid angle uj can travel unabsorbed and escape 
the galactic disk, while directions outside the solid angle 
UJ (indicated with £2 in the figure) intercept a fraction of 
the gas disk thus blocking the escape of ionizing photons. 
Note that the solid angle uj increases as the shock radius 
R{t) propagates outward. 

We can compute the solid angle uj = 27r (1 — cos 0/2) 
in terms of the opening angle (the aperture of the inner 
cones in fig. 1), which in turn is related to the vertical 
scale height h of the disk and to the radius of the ex- 
panding shell R{t) by the relation cos 6/2 = h/R{t). Thus 
we obtain the time-dependent escape fraction f{L^t) = 
UJ / 271 = 1 — h/R{t)^ where the shock radius R{t) depends 
on the AGN luminosity (see eq. 3 and text below). Since 
we are interested in the statistical, effective escape frac- 
tion, we compute the average over the duration r of the 
AGN activity 



dt 



R{t) 



(4) 



Here ts is the time (computed from the start of the AGN 
phase) at which the shock radius first encompasses the 
width of the disk (for earlier times escape fraction is null) , 
so that R{ts) = h. Using eq. (3) we express ts as a 
function of the disk dynamical time t^ as is = ts/td = 
{h/vd) / M.' Performing the integral in eq. (4) then yields 



(/(i)) = i 



ts 



1 +ln 



ts 



(5) 



where f = r /td ^ I (for interaction driven AGNs). Note 
that for shocks with M ^ I (i.e., large AGN injected en- 
ergies AE see eqs. 2 and 3) imply ts <C f, yielding 
large average escape fractions {f{L)) ^ 1, while slowly ex- 
panding shocks yield {f{L)) ^ 0. Thus eq. (5) results 
in an effective escape fraction increasing with the AGN 
luminosity L. At fixed L the escape fractions in eq. (5) 
increase with redshift since - at given - the smaller 
host galaxy masses yield larg er AE/E at higher z. 

For a ratio h/rd = 1/15 (see Narayan, & Jog ( 2QQ2[ )) de- 
termining the timescale ts for each AGN luminosity (cor- 
responding to h ^ 200 pc for a typical L galaxy), we ob- 
tain low-redshift {z o:^ 0.5) escape fractions increasing from 
(/) ^ 0.1 for M1450 = -22 to (/) ^ 0.4 at M1450 = -24, 
while at higher redshifts {z o:^ 4) the values range from 
(/) ^ 0.3 for M1450 = -22 to (/) ^ 0.8 at M1450 = -24; 
quantitative informations on the dependence of the escape 
fraction on the AGN luminosity and redshift can be found 
in fig. 2. Note that the value of the parameter h/rd = 1/15 
adopted here to compute (/) fr om eq. (5) is the same 



adopted in Menci et al. (2008), where we showed that 
the blast-wave model for AGN feedback provides column 
densities Nh (computed after averaging over all the line of 
sights schematically shown in fig. 1) depending inversely 
on the AGN luminosity. Such a behaviour is due to the 
fact faster shocks correspond to a lower fraction of still 
unperturbed gas outwards of r = Rs{t)^ and is in quanti- 
tative agreement with observations, as we showed in the 
paper mentioned above. 



3. THE AGN UV LUMINOSITY FUNCTIONS 

To predict the UV ionizing emissivity produced by our 
AGN evolution model we first evaluate the evolution of 
the AGN luminosity function as a function of redshift at 
A = 1450 A, where most of the UV rest frame data on 
AGNs are collected and where the average AGN emission 
shows the peak in the spectral energy distribution, repre- 
sentative of the overall bolometric emission. The contri- 
bution of each luminosity to the ionizing background will 
be the computed by convolving such luminosity functions 
with the escape fraction in eq. (5). 

In Fig. 2 we show the evolution of the UV luminos- 
ity function 7V(Li45o) computed in our R-SAM (binned 
in magnitude, dashed line), and the fraction A^(i^i45o) (/) 
that will contribute to the ionizing background (solid 
line), where (/) is the luminosity-dependent escape frac- 
tion given in eq. (5). The bolometric luminosity L (com- 
puted for each AGN in our R-SAM after eq. 2) is related 
to the emission L1450 at 1450 A after the bolometric cor- 



rection given in ^Marconi et al. ( 2004^ ) and Marconi et al. 
(in preparation). 

In the figure we compare our predictions with data from 
literature. Note that the data should be compared with 
the global UV luminosity function (dashed line). At red- 
shifts z < 2.5 we have compared our predictions w i th the 
luminosity function 



Richards et al 
data come from 



(|2009) 

tHe 



provided by 
and CroomeTaTn (2009). T 
combined 2dF-SDSS sample of Quasars photometrically 
selected from Sloan Digital Sky Survey, and then spec- 
troscopically confirmed using the 2DF instrument on the 
Anglo Australian Telescope. The luminosity functions are 
complete down to M1450 = —21.5 and show a gradual fiat- 
tening below M1450 = —24. The shaded region brackets 



the determinations derived in the X-ray band by La Franca 



etal. (2005); Ebreroet al. (2009); Aird et al. (2010), and 



thus accounts for systematic errors onlv. Conversions to 
1450 A have been o btained using the Marconi et al. ( 2004 ) 
relation. Note that |La Franca et al. | (|2QQ5 ); Ebrero et al 



If 

wh' 



2009) determinations account for the obscured sources 



ose fiux is reduced below the fiux limit by interven- 
ing column densities, they provide densities slighly higher 
than those of Aird et al. 2010, which are not corrected for 
this effect. 

In this redshift range, the agreement with the model 
predictions is good, if we consider that the errors shown 
are only Poissonian and that various systematics (e.g. in 
the fiux corrections) probably would widen the uncertain- 
ties involved. In any case, the luminosity functions derived 
from X-ray surveys (see shaded regions in fig 2) constitute 
upper boundaries that are fully consistent with our pre- 
dictions. 

At higher redshifts z > 3 we com pare with the rest- 
frame UV luminosity func t ions f r om [Fontanot et al. 



(2007), Siana et al. (2008) 



( 2007L iBongiorno etal. 

ana kman et al. | ( |201ip . Here the comparison with 
data at the faint end is more critical, since the corrections 
for incompleteness of the observ ations play a more rel evant 
role. For example, the data by Siana et al. (2008) have 
been derived combining IR SPITZER SWIRESata and 
optical through a standard morphological criterion and 
a color selection based on the UV dropout expected at 
these redshifts coupled with an IR color selection aimed 



Giallongo et al. 
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at removing low redshift interlopers. These authors at- 
tain a completeness level at the faintest survey limits of 
about 75% which is taken into account in their effective 
volumes. However, their completeness estimates are based 
on detailed simulations involving both templates of quasar 
spectra and a model for IGM absorption, which is crucial 
for any UV dropout selection. In their analysis they re- 
fer to the Madau (1995) model, whose analysis has been 
superseded since larger statistical infor mation about th e 
IGM absorption became available (e.g. Meiksin (2006), 
Inoue & Iwata ( [2008| ) ) , so that now the Madau model 
represents an upper limit to the average absorption. More 
recent and reliab le estimates by [Prochaska, Worseck, fc] 
O'Meara (2009) show average Lyman limit opacities by 
the IGM observed in 1800 Sloan quasars at z ~ 3.5 — 4.5 
about a factor two lower. This should allow the presence of 
a significant number of z ^ 3 quasars bluer than expected 
by the UV dropout threshold adopted by |Siana et al. | 
(2008), implying a lower completeness level of the lumi- 



nosity function di fficult to quantify . A sim ilar bias could 
be present in the "Fontanot et al. ' f2007') data which go 
two magnitudes fainter at a slightly higher average red- 
shift, z ~ 3.7. The luminosity function has been derived 
from combined standard color selection and X-ray match- 
ing in the GOODS fields. The color selection was effec- 
tive in the interval 3.5<z<5.2. In Fig. 2 we show the 
luminosity function derived in the 3.5 < z < 4 interval 
where the completeness of the survey is higher. Even in 
this survey the average IGM absorption adopted to de- 
rive the color threshold is similar to that predicted by the 
Madau model at least up to z ~ 4.5, although they apply 
some correction factor to the Lyman a forest. The impor- 
tant point in this respect is that at z > 3 the incidence 
of IGM absorption systems which dominate the Lyman 
limit opacity has a slower evolution compared to the Ly- 
man a lines (Prochaska, Worseck, & O'Meara 2009). A 
further important source of incompleteness at magnitudes 
^1450 ^ —24 is constituted by the morphological selection 
criterion adopted to reduce in the spectroscopic follow up 
the contamination from Lyman break ga l axies at the same 



redshift as described in Glikman et al. (2010, 2011) 



In the same fig. 2 we include two faint points derived 
from a recent very deep search for AGNs in the GOODS 
south field by Fiore et al. ([2012). The latter authors 
adopt different criteria starting from the standard color 
dropout technique but including faint detection in the 4 
Ms Chandra images in the 2-10 KeV X-ray band and re- 
laxing in this way any morphological restriction. They use 
spectroscopic redshifts where available or accurate photo- 
metric redshifts from the GOODS MUSIC and GOODS 



ERS multicolor catalogs (Grazian et al. 2006, 2011). In 
this way they can select faint AGN activity among the 
numerous Lyman break galaxies in a wide redshift inter- 
val 3 < z < 7 putting interesting constraints on the faint 
end of the luminosity function at z = 3 — 4,4 — 5, and 
z > 5.8 dominated by Seyfert-like objects. Although a 
fraction of the order of 20% of the detected AGNs appears 
strongly absorbed in X-ray or Compton thick, the same 
objects do not show any reddening in the optical-UV part 
of the spectrum which shows emission lines with various 
widths. Indeed subsequent spectroscopic optical follow up 
revealed even in the case of Compton thick sources, typ- 



ical highly ionized UV emission lines like CIV and NV. 
In summary, the limited spectroscopic information avail- 
able coupled with the analysis of the overall spectral en- 
ergy distributions are suggesting the absence of any sig- 
nificant link between the presence of strong absorption in 
X-ray and emission of UV photons from the central nu- 
cleus, at least in cases where the X-ray flux is measured 
for optically-UV selected AGNs. In this respect the sample 
is more complete respect to the optical surveys a nd this is 
reflected in the higher densities attained by the [Fiore et 
al. I (20121) luminosity functions. Nevertheless the sample 
could be affected by several uncertainties and significant 
incompleteness, especially at the faintest limits where the 
objects have been selected from very deep HST NIR im- 
ages at H ^ 27 and measured at the faintest X-ray fluxes 
of Fx ~ 2 X 10~^^ ergs s~^ cm~^. At the lowest X-ray 
luminosities, objects with relatively high X/Optical ratio, 
and consequently with H > 27, are missed in their survey. 
Most of the redshifts have been evaluated photometrically 
and this could introduce low redshift interlopers in the 
sample. On the other hand, even in this case the color 
selection based on the standard UV dropout technique 
adopts a predicted IGM opacity hi g her th an measured by 



Prochaska, Worseck, & O'Meara (2009) although lower 
than expected in the Madau model, allowing for some in- 
completeness. As a consequence uncertainties in the Fiore 
et al. LF data due to systematic errors could be larger than 
shown in the figure based on number statistics. Keeping 
in mind these limitations we have adapted their LF data 
after having recast the X-ray luminosities in UV (1450 A) 



(2004) 



magnitudes following the standard |Marconi et al. 
corrections. 

At redshifts z > A the above observational biases and 
possible incompleteness could explain the model over pre- 
diction appearing at fainter magnitudes M1450 > — 24 re- 
spect to the present dataset. The discrepancy is a factor 
2 — 3 in the magnitude interval — 24 < M1450 < —22 up 
to z ~ 6 and increases to fainter magnitudes. Since the 
estimate of the observational incompleteness is difficult to 
assess without the data input from deep X-ray surveys, we 
shall rely on our model to get an insight on the contribu- 
tion of AGNs to the ionizing background. To this aim, we 
shall couple our predicted luminosity functions with the 
physical description of the escape fraction based on the 
blast- wave model. 

4. THE AGNS LYMAN CONTINUUM EMISSIVITY AND 
HYDROGEN PHOTOIONIZATION RATE 

We now proceed to compute the contribution of our pre- 
dicted luminosity functions to the AGN emissivity at the 
Lyman limit. This is given by 



{f{L))N{L)L^dL 



(6) 



where the {f{L)) is the fraction of AGNs of a given lumi- 
nosity which shows Lyman Continuum emission given by 
eq. (5), and E/Egi2 is the (luminosity dependent) bolo- 
metric correction at 912 A. In the integral we have consid- 
ered all the AGNs produced in our MonteCarlo code down 
to very faint magnitudes, however the total emissivity at 
a given redshift is independent of the faintest luminosity 
cutoff. This is clearly shown in fig. 3, which shows the 
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contribution to the total emissivity from AGNs with dif- 
ferent bolometric luminosity L at different redshifts. Since 
the bolometric correction from the UV M1450 magnitude 
is not linear, we have indicated on top of the first panel 
the UV magnitudes at 1450 A corresponding to the peak 
bolometric luminosities. Note that while at 2; ~ 3 the 
main contribution comes from AGNs in the relatively wide 
magnitude interval — 24 < M1450 < —18.5, at higher red- 
shifts the main contribution comes from brighter AGNs 
with — 24 < M1450 < —22. In this magnitude interval the 
discrepancy between data and model prediction is within 
a factor 2 — 3, which is of the same order as the differ- 
ences among surveys where AGNs are selected with differ- 
ent methods in various bands. 

In this respect it is just in this magnitude interval that 
a more careful analysis should be performed to investi- 
gate whether data and model prediction can be reconciled. 
Unfortunately this magnitude interval is sparsely sampled 
by present surveys; more statistics with much more con- 
trolled systematics should be produced in the future. In 
particular, multiwavelength surveys of medium depth (e.g. 
/ ^ 22) in an adequate sky area should be planned for 
spectroscopic and X-ray follow up. In the meanwhile we 
explore the ability of our linked AGN/galaxy model in 
providing sufficient UV photons to keep the IGM ionized 
up to z ~ 6. In fig. 4 and Table 1 the predicted 
ionizing UV emissivity 624 in units of 10^^ erg s~^ 
Hz~^ Mpc~^ is shown as a function of redshift. 
The sharp increase up to z ~ 2 is due to the lu- 
minosity evolution of the AGN luminosity func- 
tion. At z > 4.5 the emissivity gradually declines 
up to z ~ 8 and then drops quickly. The rather 
flat behaviour at high redshifts is due to the fre- 
quent galaxy merging rate which drives the AGN 
activity coupled with the short gas cooling time. 
Both guarantee a high AGN duty cycle and conse- 
quently an average AGN emissivity weakly evolv- 
ing in time. 

The photoionization rate per hydrogen atom r_i2 in 
units of 10~^^ s~^ is the standard quantity used to eval- 
uate the ionization status of the IGM. We have com- 



puted this quantity following |Madau, Haardt, fc Rees 
(T999| and |Schirbe'r fc Bullock | ( |2QQ3p 



r_i2(z) c^O.46 



624 (^) / A/ \ 

3+ I auv I \bOMpcJ 



4.5 



3-(2.5+7) 



(7) 

where auv = —1.76 is the average UV spectral index of 
AGNs derive d from the average SED we adopt (Marconi' 
et al. 2004). The power-law index —(2.5 + 7) describes 
the decrease in redshift of the mean free path (mfp) Al of 
ionizing photons in the IGM due to the increase in red- 
shift of the Lyman Limit absorption systems with slope 7. 
We adopt 7 = 1.94 and the normalization of t he mfp to 50 
Mpc at z = 3.5 as in Songaila fc Cowie | ( [2010D . At z > 4.2 
we adopt a steeper redshift evo lutionwith 7 = 5 to fit the 
mfp data point at z :^ 5.8 (see Songaila fc Cowie (2010) 



for details). This equation represents a first order approx- 
imation where only ionizing sources within one absorption 
length (where t/gm = 1) contribute to the ionizing back- 
ground. The resulting curve from our MonteCarlo model 
is shown in Fig. 5. The dotted one represents the intrin- 
sic production of UV ionizing photons by the overall AGN 



population. The ratio between the two curves depends on 
the ratio between the Lyman limit optically thin AGNs 
and the overall population. 

Data points with different symbols are also shown for 
comparison. They represent different estimates derived 
in a model-dependent way from the IGM statistics. In- 
deed most values come from methods of estimating the 
UV background; and so F based on the mean flux decre- 
ment in the Lyman a forest of QSO spectra in combination 
with numerical simulations. The UV background value is 
adjusted until the mean flux in si mulated Lya forest spec- 
tra is equal to that in real data (Faucher-Giguere et al. 

2008; Boltor TeFlir" 2005; Beck er, Ranch, fc SargenT 
,2007, Fan et al. |2006[ tWyithe fc Bolton |201ip . An alter- 
native method is based on the proximity effect which has 
recently been pro posed at the highest redshifts z = 5 — 6 by 
Calve rley et al. ||201 1). However both methods are sub- 
ject to uncertainties and systematics which are responsible 
for the scatter among the data points reported in Fig. 5. 
The data shown in the figure derived from flux decrement 
analysis in the IGM Lyman forest of QSO spectra, have 
been scaled to the sam e IGM temperature-density relation 
as in Calverley et al. (2011). 

Indeed, converting the mean flux decrement into an ion- 
ization rate depends on modelling the gas density and 
temperature distribution in the IGM at very low densi- 



ties (Miralda-Escude, Haehnelt, fc Rees [2000 ) in various 
cosmological scenarios. On the other hand, estimates from 
the proximity effect are sensistive to the assumption con- 
cerning the presence of overdensities of matter close to the 
measured quasar and to the assumption that gas temper- 
ature within t he proximity region is similar to that in the 
general IGM j Calverley et"al|[2QTT ). 

The continuous curve represents the hydrogen ionization 
rate as a function of redshift predicted by our model, where 
the oscillating behaviour is related to the stochastic nature 
of the AGN trigger, and rendered by our MonteCarlo code. 
The model is able to reach an almost constant value of the 
order of r_i2 — 1 consistent with the data points, in a 
broad redshift interval from z ^ 2 up to z ^ 4. A decline 
by a factor 4 is then appearing at higher redshifts up to 
z ^ 6 mainly due to the significant decrease of the mean 
free path of ionizing photons. 

The gradual decline at high redshifts is a strong intrinsic 
prediction of our model based on a linked AGN/Galaxy 
evolution and is in line with recent data analyses a t 
z - 5-6 dCalverley et al. |20lH[Wyithe fc Bolton |201lD . 
As we have shown in the previous section, the model over- 
estimates the number of relatively (M1450 ~ —23) faint 
AGNs at z > 4 by a factor ^2 — 3; on the other hand we 
argued the possible presence of significant incompleteness 
in the scanty present faint dataset. 

Since the predicted emissivity by hard spectrum sources 
like AGNs keeps high up to z ~ 6 it is interesting to eval- 
uate the impact of this evolution on the ionization his- 
tory of Hell in the IGM. Indeed the higher Hell ionization 
threshold respect to hydrogen and its smaller photoioniza- 
tion cross-section imply a complete Hell reionization only 
when numerous hard spectrum sources like AGNs are dis- 
tributed within the IGM. There is some consensus about 
a late Hell reionization which could complete at z ^ 3 
as suggested by observations of Hell Lya Gunn-Peterson 
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trough and patchy absorp tion in the Hell Lya forest at 



3 (Shuh et al. (2010) and reference therein). Other 



indirect evidence comes e.g. from an apparent increase 
of the average IGM t emperature at z ~ 3 measured from 



metal line ratios (e.g. Songaila JT998)). 

Following early analysis by 'Ma dau, Haardt, fc Re es 



([1999) and more recently by Haardt & Madau I (|2012|) we 
nave computed the volume filling factor of Helll in the 
IGM as a function of redshift. We have also considered 
the case where Lyman limit absorption by optically thick 
clouds (LLSs) limits the AGN io nization, as proposed by 
Bolton, Oh, & Furlanetto (2009|). In the latter case each 
AGN is able to ionize only a limited region due to the fact 
that, at variance with hydrogen ionization, the distance 
between Hell Lyman limit optically thick clouds is smaller 
than the typical AGN separation at z > 3. Thus Hell- 
ionizing photons will have intersected several helium LLSs 
being subject to absorption. The equation describing the 
ev olution of the Helll volume filling factor Q is, according 

RQ - fiQ 



Bolton, Oh, & Furlanetto (2009): Q 



to_ ^ 

where h is the Hell ionization rate per helium atom and 
R{t) is the average recombination rate which includes an 
effective clumping factor. The last term describes absorp- 
tion by LLSs clumps. In Fig. 6 we show the evolution of 
the volume filling factor assuming only recombination with 
moderate clumping factors. We also show the evolution 
including LLSs absorption as a further photon sink as in 
Bolton, Oh, fc Furlanet'tQl ( |2009D . The Q evolution adopt- 
ing the 'Haardt fc Ma dau"] ( |201 2) emissivity is also shown 
for comparison. It appears that the Hell reionization can 
be completed down to z ~ 3 assuming modest clumping 
especially when LLSs allow each AGN to ionize only a 
limited volume. The early hydrogen reionization is shown 
for comparison only in the recombination case since the 
mfp of hydrogen LLSs is larger than typical AGN sepa- 
ration and the photon sink term can be neglected. Of 
course the knowledge of the detailed Hell reionization his- 
tory requires 3D hydrodynamical models. The adopted 
simplified description is only suggesting that some degree 
of inhomogeneity could help to reduce the tension between 
the possible presence of a high level of ionizing emissivity 
by AGNs at z ~ 6 with an Hell reionization ending at 
z ^ 3. 



5. SUMMARY AND CONCLUSIONS 

A faint AGN population is emerging from recent multi- 
wavelength surveys, in particular from those including X- 
ray detections. These are changing our perception about 
the AGN contribution to the ionizing UV background, al- 
though the selection of faint AGNs at the highest redshifts 
z b — 6 becomes difficult with the present instrumenta- 
tion. Given the current observational limits, insight into 
the origin of ionizing ffux at high redshifts can be gained 
from theoretical modelling. We have adopted our semian- 
alytic model R-SAM which successfully links the evolution 
of the galaxy population to that of the AGN population 
through the computation of the growth of supermassive 
black holes in the nuclei of galaxies. The model was able 
to reproduce several observables such as luminosity and 
mass functions of both galaxies and AGNs. 

In this context we first concentrated on the compari- 
son between the predicted and observed LFs at the UV 



rest-frame wavelength 1450 A. The agreement is satisfac- 
tory at low and intermediate redshifts up to z ~ 3. At 
higher redshifts the predicted LFs tend to overestimate 
the observed data by a factor 2 — 3 at the fainter mag- 
nitudes Mi45o ~ —22. We have discussed the possible 
reasons for the discrepancy. From the observational point 
of view it is clear that surveys of optically/UV selected 
AGNs which also rely on X-ray detection as a discrimi- 
nant in the candidate selection process appear more com- 
plete on approaching the predicted values by the model. 
Indeed without the X-ray discriminant, the point-like mor- 
phological criterion is needed to avoid confusion with the 
more numerous Lyman break galaxy population. But this 
introduces significant incompleteness at faint magnitudes. 
Moreover several uncertainties in the selection process like 
e.g. the lack of spectroscopic information for the faintest 
AGN candidates, make the estimates of the faint end of 
the AGN LF appreciably uncertain at very high redshifts. 

Since faint AGN surveys could be affected by significant 
bias we have relied on our model to get an estimate of the 
hydrogen photoionization rate produced by the AGNs. To 
this aim, we have coupled our predicted luminosity func- 
tions with the physical description of the escape fraction 
based on the blast-wave model already introduced in our 
R-SAM to describe the expansion of the blast-wave pro- 
duced in the interstellar medium of the hosting galaxy by 
the AGN-driven outfiows. The blast wave clears the way 
for the UV photons emitted by the nuclear AGN which 
are free to escape outside the galaxy and ionize the inter- 
galactic medium. To evaluate whether AGNs can ionize 
the intergalactic medium we have first computed the AGN 
emissivity at the Lyman limit as a function of luminosity 
and redshift. Interestingly, the predicted AGN emissivity 
has a definite maximum at about M1450 ~ — 23 which is 
almost independent of redshift. Fainter sources give a neg- 
ligible contribution because of the fiattening of the shape 
of the luminosity function and because the escape fraction 
of UV photons decreases on average in fainter AGNs and 
reaches (/) ^ 0.1 — 0.3 at M1450 ~ —22 at intermediate 
and high redshifts, respectively. Note that for AGNs of 
a given luminosity the escape fraction increases with red- 
shift following the increase of the merging rate which is the 
main triggering of the QSO activity in our model. Since 
at z > 3 AGNs with M1450 — 22 show a duty cycle of 
about 10% and so represent a bout 10% of the host galaxy 
population ( Fiore et al. [2012 ) then the bulk of the Lyman 
Continuum emitters hiding faint AGN activity represents 
about 1% of the UV-dropout galaxy population similar 
to recent statistics on large Lyman break galaxy samples 
( jVanzella et al. ,2010^ . The resulting photoionization rate 
evolves strongly from z = to z = 2, and then remains 
nearly constant to a value r_i2 ~ 1 up to z ~ 4.5. A 
gradual decline is apparent at higher redshifts, consistent 
with the trend derived from most recent analyses of the 
Lyman forest in QSO spectra. We have also shown that 
a high level of AGN ionizing emissivity up to z ~ 6 could 
still be consistent with a more gradual Hell reionization 
ending at z ~ 3, especially if absorption by Hell Lyman 
limit IGM clouds is taken into account. 

In summary, the main prediction of the model is that in 
a scenario where the AGN activity is linked to the dynam- 
ical galaxy evolution triggered by interactions and merg- 
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ing, only a mild high redshift evolution of the hydrogen 
photoionization rate is expected in the redshift interval 
z ~ 2 — 5. New deep multiwavelength surveys with ac- 
curate selection procedures are needed to build reliable 
AGN volume densities down to M1450 ~ —23 in the rele- 
vant redshift interval. These values are critical to under- 
stand whether AGNs can play a key role in reionizing the 
high redshift IGM. Our model suggests that it is time to 
reconsider the AGNs as the main population driving the 



ionization history of the Universe over the whole known 
redshift interval (up to z ^ 6) where the cosmic evolution 
of structures is observed. 

We thank the referee for useful comments which have 
significantly improved our paper. We are grateful to A. 
Cavaliere for critical reading of the manuscript. We also 
thank F. Haardt for useful discussion on the Helium reion- 
ization history. 
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Fig. 1. - A schematic representation of the geometrical effects determining the escape fraction in the blast wave model of AGN feedback. The 
shock radius Rs{t) expands outwards, compressing the swept gas into a thin shell (represented in darker colour) with width — A) (see 

, Cavaliere, Menci | ( [2b05| for a computation of A), and leaving a cavity inside. Such a density distribution p(r) is plotted in detail at 
the bottom. The solid angles uj subtend the intercept of the expanding blast wave with the galactic disk, and are represented as grey (yellow 
in the electronic paper) cones. Ionizing photons can travel unabsorbed only along the directions contained within a;, while the directions 
external to the cones (indicated as -^2) intercept a fraction of the galactic gas, thus inhibiting the escape of ionizing photons. 
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Fig. 2. - The evolution of the UV luminosity function obtained from our SAM model (lines) is compared with different observational data. 
Dashed lines refer to the full intrinsic emission at 1450 A, while solid lines show the fraction of objects that will contribute to the ionizing 
background and suppress the luminosity function by the luminosity-dependent escape fraction given in eq. (5). The redshift intervals and the 
data we compare with are shown in the different panels and discussed in the text. 
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Fig. 3. - The integrated emissivity of the AGN population (computed from eq. 6) contributed by AGNs with different bolometric luminosity 
(on the X-axis) in different redshift bins. The histograms represent the differential contribution, while the cumulative contribution is represented 
by the dotted line. The UV magnitudes at 1450 A corresponding to the peak bolometric luminosities are also shown on the scale at the top of 
the top-left panel. 
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Fig. 4. - Lyman continuum emissivity produced by the AGN population in the R-SAM model as a function of redshift. The curve has been 
smoothed respect to the original MonteCarlo values. 
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Fig. 5. - Photoionization rate per hydrogen atom produced by the AGN population in the R-SAM model as a func tion of redshift is compared 
with data from different authors scaled to the same IGM temperature-density relation as in'Calverl ey et al~| ( |2011| >. The solid line is the model 
prediction after eq. 6 and 7, while the dotted one represents the intrinsic production of UV ionizing photons by the overall AGN population 
(computed from eq. 6 and 7, but assuming an escape fraction (/} = 1 in eq. 6). References for the different data points are shown on top of 
the figure. 
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Fig. 6. - Volume filling factor Q of Helll and HII as a function of redshift. Dashed (red) curve shows Q(Helll) with He recombination and 
clumping factor C = 10; long dashed (viole t) curve shows Q(Helll) with C = 7; continuous (black) thick curve also includes photon sink by 
LLSs and C = 3] [Haardt Madau~| ( |2012| ) is shown as dotted (blue) curve for comparison; thin continuous (black) curve is the HII volume 
filling factor evolution assuming C = \JT. 



